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Abstract 

The complete matrix elements for e^e~ bbW^H^ and e~^e~ bbH^H^ 
are computed at tree-level within the Minimal Supersymmetric Standard 
Model. Rates of interest to phenomenological analyses at the Next Linear 
Collider are given. In particular, we analyse: (i) ti pair production followed 
by the top decays into bbW^H^ and bbH~^H~ final states; (ii) H^H~ signals 
in which hW^ bbW^ and ^ tb ^ bbW^. Top and Higgs 

finite width effects are included as well as all those due to the irreducible 
backgrounds. 
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1. Introduction 



In the present paper we intend to conclude a long-term project of computing effects 
due to the finite width of the unstable particles as well as to the non-resonant diagrams 
entering in electroweak (EW) processes of the type 

e+e" ^ bbX+X-, (1) 

where = W^, H^, within the Standard Model (SAi) and the Minimal Supersym- 
metric Standard Model {A4SSA4) at Next Linear Collider (NLC) energies [|l|. The 
case X~^X~ = W^W~ within the SM was considered in Refs. [0, g (see also Ref. IQ) 
whereas the equivalent process in the M.SSM. was discussed in Ref. 00. 

By the time that an NLC machine will be operating, will already be known pre- 
cisely and either will have been measured or a large part of the M.SSM. parameter 
space [Ma-, tan /5) defining the values of Mniggs will have already been explored^, thanks 
to the experiments at the Tevatron , LEP2 and LHC |10| . (For an illustration of 
the discovery potential of such machines, we refer the reader to 0, ||, |lOl-) Anyhow, 



whichever the status of the measurements of m^, Mu-, M/i, Ma and/or M}i± will be, we 
would like to stress that detailed studies of both top quark and Higgs boson properties 
(other than the mass, e.g., the width, the branching ratios, etc.) will more easily be 
carried out at an e^e~ linear accelerator, operating in the energy range Ivnt ~ a/s ~ 1 
TeV. At such a machine, the absence of the huge QCD background typical of hadron 
colliders combined with the large luminosity and 'effective' centre-of-mass (CM) energy 
available allows one to cover the whole of the Higgs mass range and to exploit all possi- 
ble Higgs and top decay channels in order to perform very high precision measurements 
of their parameters. 

Top and Higgs signatures naturally enter in processes of the type ([l|). In fact, in 
the SM.^ the process e^e~ bbW^W~ is a signature of both top quark production 
in ti pairs and Higgs boson bremsstrahlung in the Z(f) channel. On the one hand, top 
pairs produced via the e^e~ ti decay through ti — > bbW^W~ whereas, on the other 
hand, the channel Zip bbW^W~ will be one of the cleanest ways of detecting a 
heavy Higgs, thanks to the high performances expected from the vertex detectors in 
triggering 6-quarks pT|p|. Within the AiSSJ^, the reaction e~^e~ bbW~^W~ allows 
one to study not only the Supersymmetric counterparts of the two mentioned SAi 
processes^, but also AH production followed hj A ^ bb and H W^W~ as well as 
hW+W- in which h ^ bb ^. 

It is the purpose of this note to consider the last two possibilities left within the 
MSSM, that is 

e+e" ^ bbW^H"^, (2) 
e+e" ^ bbH+H-. (3) 

^For the case of the irreducible background from QCD in + 4jet events, see Ref. ||]. 

■^Here and in the following, (p indicates the Higgs boson of the SAi whereas the subscript 'Higgs' 
refers collectively to the H, h, A, scalars of the AiSSAi. 

"^The branching ratio (BR) of the channel Z — > 65 is five times larger than that into /^^/i^ or e+e^ 
and it is free from backgrounds due to decays. 

^In case of Higgs production, the heaviest of the A4SSA4 neutral Higgses (i.e., H) replacing the 
SM. Higgs scalar (i.e., 0). 



As it was done in the previous publications 0^0; our aim is twofold. 

• Firstly, to establish the importance of top finite width effects and of those due to 
the non-resonant background in case of ti bbW^H^ and tt bbH^H^ decays. 

• Secondly, to study Higgs signatures, such as e+e^ — > H^H production followed 
by the decay of one of the charged Higgses into: 

1. h*W^* bbW^* and 

2. t*b bbW^ (if mt ^ Mh± - rrib), 

where the asterisk indicates possible off-shell channels. Note that these two signa- 
tures are produced via process @. We further notice that the top-antitop decay 
channels 

3. ti bbW^H^ and tt bbH+R- 
represent themselves a source of iJ^ signals. 

Although all the topologies yielding bbW^H^ and bbH~^H~ final states produce 
Higgs bosons (both charged an neutral), we will consider the corresponding diagrams 
as (irreducible) 'background' to the two 'signals' e~^e~ tt and e~^e~ — >• H~^H~, as 
the former represent Higgs production mechanisms that are largely suppressed with 
respect to the latter and also to others exploited in neutral Higgs boson searches (for 
a review of these see, e.g., Ref. [|12|)|^. Furthermore, it is worth reminding the reader 
that, at least over the areas of the AiSSAi parameter space where the t — > bH~^ and 
—>■ hW^ — > bbW~^ decay modes are both important (see later on), the processes 
e'^e~ ti and e~^e~ —>■ H^H~ cannot be unambiguously separated and studied inde- 
pendently. Therefore, in such (MA,tan/3) regions, any of the two reactions constitutes 
an irreducible background to the other. The different topologies contributing at leading 
order to processes (0)-(^ are given in Figs, la and lb, respectively. In the first case 
they lead to 51 different Feynman diagrams whereas for the second case the number is 
60. 

For future reference, we also present in Fig. Ic the M.SSM. BRs of top quarks 
and charged Higgs bosons in the Mh± range and for the values of tan/5 that will be 
used in our analysis. Note that such rates dictate the decay phenomenology of the top 
(diagrams 2 in Fig. la and 3 in Fig. lb) and Higgs (diagrams 9 and 12 in Fig. la) signals 
that will be studied later on. As in our forthcoming discussions we will interchange the 
role of Mh± and Ma-, we have plotted both these quantities as mass scales in Fig. Ic. 

The plan of the paper is as follows. In the next Section we describe the calculations 
we have performed. In Section 3 we present our results. Our conclusions are given in 
Section 4. 

^We notice that such diagrams arc nonetheless interesting on their own, e.g., as a way of testing 
the gauge structure of Higgs sector of the as they involve many of the fundamental tree-level 

vertices of the theory. This is however beyond the scope of this study. 



2. Matrix Elements 



To compute the Matrix Elements (MEs) of processes (@)-® we have used hehcity am 
phtude techniques. A first code has been produced by using the method of Ref. 



This has then been checked against a second one buih up by resorting to the subrou- 
tines contained in the package HELAS |]14|. In particular, in order to implement the 



diagrams of Figs, la-b, we have made use of the subroutines involving scalars with 
the vertices and propagators appropriately arranged to reproduce the interactions and 
wave-functions of the five Higgses in the unitary gauge of the J^SSA4. 

The two FORTRAN codes have been programmed by the two authors independently 
and they agree at amplitude squared level within 12 significative digits in REAL*8 
precision. Moreover, all the amplitudes produced with the first method have been 
tested for gauge invariance by computing the two MEs in other two gauges as well (i.e., 
the Feynman and the Landau ones). In order to deal with the complicated resonant 
structure of processes (@)-(|l), we have adopted the technique of sphtting the two MEs 
in a sum of non-gauge-invariant pieces, each of these implementing a different resonant 
structure, and of integrating them separately with the appropriate mapping of the 
phase space variables. Such a procedure has already been documented in detail for the 
case of the e^e^ — > bbW~^W~ process in the SM. in Ref. 0, so we only sketch out here 
the method. 

Firstly, one isolates the diagrams with similar resonant structure by grouping these 
together in 'subamplitudes'. For example, one can recognise in process (|^) the reso- 
nances: ti {bW+){bH-), t bW+, t bH-, H, h,A,Z ^ bb, H+R- {tb)H- 
{bbW+)H- and H+R- {hW+)H- {bbW+)H-. In case of reaction (D one 
has the following Breit-Wigner peaks: ti {bH^){bH^), t — bH^, i bH^ and 
H, h,A,Z^ bb. Secondly, one defines the mentioned non-gauge-invariant components 
of the amplitudes squared, by appropriately combining the subamplitudes. For exam- 
ple, in both processes (@)-(|D we have simply taken the square of the various resonant 
subamplitudes (yielding 9 and 7 terms, respectively) whereas we have combined their 
interferences altogether in the last term, apart from those between the two single-top 
and the only double-top channels, which were added to the t- and t-diagrams in both 
reactions (||)-(H). Thirdly, one maps the phase spaces around the resonances. For 
example, for our analysis this was made by performing the change of variable 

Q'-M' = MT tan 9, — . dQ' = (Q^ - M^r^ ^^ 

(where Q, M, T stand for the virtuality, the mass and the width of the resonance), 
which gives integrands smoothly dependent on 6. Fourthly, the various amplitude 
squared terms are integrated separately (we have used VEGAS to this purpose) 
and added up in the end (to recover gauge invariance). We however notice that, if 
taken separately, these terms provide a useful way of looking inside the processes and 
distinguishing between the different fundamental interactions. 

Concerning the values adopted for the SA4 and AiSSA4 parameters and the rela- 
tions implemented for the AiSSAi Higgs masses and couplings we refer to Refs. 
Furthermore, as in the previous publications 0-0, we have chosen here, as represen- 
tative for tan/3, the two extreme values 1.5 and 30. Ma varies between 50 and 220 



GeV, as higher values would be uninteresting in the present context (the corresponding 
rates for both the t bH^ decays and the H'^H~ pair production become heavily sup- 
pressed). Following the latest measurements at FNAL |16], the top mass has been set 
equal to 174 GeV. Therefore, we use as indicative of the top-antitop threshold energy 
the value -y/i = 350 GeV, along with the final stage energy of 500 GeV. 

Our calculations have been carried out at leading order only. We have not included 
beamsstrahlung and bremsstrahlung effects either. We are however confident that such 
approximations will not spoil our conclusions. Indeed, it has been proved that this is 
true for the process e~^e~ — > hhW^W" p|, ^ and we expect the same in the present 
context. 

A final consideration is now in order. Following the approach of Refs. pl-|5|, we 
have generally kept the final state bosons and H"^ on-shell in our computations. 
However, we made an exception in one case: when computing the squared of graph 
12 of Fig. la. This represents the Higgs signal e^e~ H^H~ in which one of the 
two scalars decays into hW^ pairs, followed by /i — > hh. That is, when only the light 
neutral Higgs is allowed to be off-shell. However, this channel is relevant in the case 
of off-shell bosons as well. Therefore, in order to keep this into account, when 
calculating the contribution of the mentioned diagram to the total cross section of 
process (|^) we have attached a fermion-antifermion decay current to the boson and 
computed the amplitude squared of a 2— >5 process. Indeed, proceeding this way, the 
rates that we calculated from the integration of this diagram over the appropriate phase 
space reproduce rather accurately those obtained by multiplying the cross section for 
e+e" H+H~ times the off-shell BR[H^ hW^ where j represents 

a jet (see next Section). 



3. Results 

In carrying out our analysis we closely follow Refs. 0, ||, ^ and the studies reported 
in the Sections 'Higgs Particles' and 'Top Quark Physics' of Ref. |T]]. Furthermore, 
throughout this paper we assume that the mass scale of the Supersymmetric partners 
of ordinary matter is beyond the energy reach of the NLC, so that they cannot be 
produced at such a colhder. In particular, we neglect considering the Supersymmetric 
decay t — > ixi, where i refers to the stop and x? is the lightest neutralino, as well as 
other possible MSSM modef] 0. 

3.1 Event selection 

In our selection strategy we assume a six-particle semi-leptonic(hadronic) tagging, for 
both reactions (^-(^: that is, 6's are assumed to hadronise whereas one of the charged 
bosons decays into jet-jet pairs and the other into rur- In this way, the selected final 
state is always made up by a high energy and isolated r accompanied by a multi- 
hadronic system and appreciable missing four-momentum. 



^For a recent review of the latter, see Ref. p8{ . 



This means that in events of the type (|^) one can in principle pick up a tau-lepton 
either from a or a H"^. However, it should be noticed that the four- momentum 
of the tagged boson is eventually reconstructed by summing that of the decay lepton 
to the missing one (the latter being computed using the kinematic constraints |I^|)p|. 



Therefore, one can always ask that the invariant mass of the lepton-neutrino system 



is larger than \J M^± + M\ f» 100 GeV to remove — > ru decays. Ma ~ 60 GeV 
being the lower experimental limit on the mass of the pseudoscalar neutral Higgs of the 
AiSSAi [0. This way, we explicitly require only hadronic decays in process (0). 
This has a clear advantage: in fact, it allows one to easily reconstruct resonant invariant 
mass spectra (i.e., Mf,iY±^^^ for the top and Mf^l^r±_^^,ljj for the charged Higgs). In 
addition, the hadronic channels represent the largest fraction of the total decay rate 
(around 70%). We also believe the selection of the channel if^ — rz/^ in reaction (0) not 
restrictive in the present context, for two reasons. On the one hand, the tau- neutrino 
channel is always sizeable up to Mh± ~ 170 — 180 GeV (i.e., the tb decay threshold) 
and for all tan/5's. Jl^. Therefore, when tagging semi-leptonically(hadronically) over 
this Mh± range one naturally retains a large part of the original top (diagram 2 in 
Fig. la) and Higgs (diagram 12 in Fig. la) events. Indeed, it should be noticed that 
only charged Higgses with mass up to 161 GeV or so (corresponding to Ma ~ 140 
GeV) can be produced with detectable rate at ^/s = 350 GeV via e~^e~ H^H~ . On 
the other hand, when Mh± ~ 180 GeV (thus when ^/s = 500 GeV), top quarks mainly 
decay via the bW^ channel so that e~^e~ tt events are in practise accounted for 
by the e"'"e~ — > bbW^W~ channel alone p|, In addition, we remark that double t- 
and/or /i-decays H^H^ bbti,tbhW^,hhW~^W^ bbbbW^W^ of charged Higgses 
produced via e^e' H^H^ would yield a final signature with numerous hadronic 
tracks (e.g., up to 8 jets can be produced, if W^W~ — >• jjjj !) with the consequent 
problems due to the combinatorics, so that to tag one of the Higgses via the lepton- 
neutrino channel might be a better approach^ For example, for large tan (3 values, the 
H"^ tVt- channel has a BR that is always larger than 0.35 [|1^ for the values of 
masses that can be produced at ^/s = 500 GeV (say, up to Mh± ^ 235 GeV, that 
is, for Ma ~ 220 GeV). In contrast, if tan/3 is small, the heavy mass region Mh± ~ rrit 
is clearly sacrificed by requiring one of the if^'s to decay into rz/ pairs Jl^, so that 



in such a case the mentioned multi-hadronic final states will necessarily need to be 
exploited in experimental analyses. 



e+e- bbH+H- 

In the case of reaction (^ one is only interested in top decays (diagram 3 of Fig. lb), 
since events cannot be produced at large rate by any of the other diagrams. There- 

^This is exactly true only when the energy of the colliding electron/positron beams is known (i.e., 
no Linac energy spread, bremsstrahlung and beamsstrahlung) , when the detector is ideal (i.e., infinite 
resolution and complete geometrical coverage) and when no other neutrino from cascade decays is 
present in the interaction. In a more realistic case, corrections have to be made to account for all the 
above effects, which in fact do not spoil the effectiveness of the method [ p9[ . 

^Alternatively, such multi-hadronic channels would need a proper simulation (including hadronisa- 
tion, detector effects, etc.) which is clearly beyond the intention of this study. 



fore, in this case, the relevant charged Higgs masses are those for which Mh± < nit— nib 
and strictly below the top-bottom threshold the rvr channel is the dominant one at 
all tan/3's [|1^. In these conditions then, a purely leptonic decay channel might seem 
more appropriate: i.e., ti hhH^H^ —>■ bbT~^ t~ Uri^r ■ However, we do not consider 
here this possibility. Firstly, because it would prevent one from reconstructing invariant 
mass top/Higgs spectra and, secondly, because over the experimentally allowed tan/3 
range (say, between 1 and 40) the t —>■ hH^ decays are always suppressed compared to 



the t —> bW ones |21|, so that the tt —>■ bbW channel represents a better option 



anyway to detect charged Higgs scalars 



Our results are presented in Tabs. I-II and Figs. 2-80. 
3.2 Top signals 

As the techniques to measure the top parameters differ greatly depending on whether 
the collider is running around the 2mt threshold or far above that, we will discuss 
separately the main features of processes at the two CM energies considered. 

= 350 GeV 



Since one of the means of measuring mt at an NLC is a threshold scan ||T9[, Figs. 2a 
and 2b show the total cross sections in an interval around the top-antitop threshold at 
2mt for the case of bbW^H"^ and bbH~^H~ final states, respectively, for a selection of 
pseudoscalar Higgs masses (for larger values of Ma top decays into charged Higgs are 
well below detection level): in the case of the e^e~ —>■ tt diagram alone as well as in 
presence of all graphs. Whereas in the case of two charged Higgses (see Fig. 2b) the 
irreducible background is almost negligible (only a ten percent effect at the most can 
be noted very below threshold, at small Ma's), when only one of the top quarks decays 
into bH^ pairs its effect can be drastically visible (see Fig. 2a). This is always true 
at small tan j3, whereas for large tan /3's this is the case only at smaller values of Ma 
(for Ma = 140 GeV and tan (3 = 30. the ti and the complete bbW^H^ rates virtually 
coincide). To appreciate quantitatively the various effects entering in the curves of 
Fig. 2a-b we have displayed in three different columns (see Tabs. I and II) the cross 
sections due to: (i) the e~^e~ — > ti diagram in Narrow Width Approximation (NWA, 
r 0)0; (ii) the e~^e' —>■ ti diagram with a finite width (F = F^); (iii) all the diagrams 
(again with F = F^) for both processes (0)-(§) and the two tan/3 values. This allows 
one to isolate the effect of the finite width F^ and of the irreducible background. We 
also report rates for three values of Mh± so that one can notice the heavy suppression 
of the total cross sections, which sets up as the charged Higgs mass increases and which 

^°Please note that, when reporting the rates corresponding to process (H) in the forthcoming Figures 
and Tables, we have considered only one of the two possible charge combinations of the and 
bosons. The cross sections for e+e" hhW'^H~ and e^e~ — > hhH^W~ are in fact identical because 
of charge conjugation. Therefore, the total production rates correspond to twice the values given in 
Tabs. I-II and Figs. 2-8. We have decided to do so in order to be able to plot in the same figures 
the spectra of both processes the rates of the latter being generally smaller than those of the 

former. 

"See Refs. [i, i for details. 



is due to the decreasing branching ratio of the top into bH^ pairs. 

Thanks to our procedure of computing the 'total cross sections' (i.e., by summing 
different 'sub-cross sections' with a specific resonant structure) we have been able to 
identify which of the diagrams in Figs, la-b are responsible for the larger background 
contributions entering in Figs. 2a. 

For example, for tan/3 = 1.5, at Mj^ = 60 GeV, h{A) bb resonant graphs (not 
including graph 12 in Fig. la) yield a contribution of ~ 0.16(0.41) fb, at y/s = 340 
GeV. Of the two Higgs signals only — »• hW^ —>■ bbW^ is relevant, producing a 
rate of ^ 8.68 fb, whereas tb ^ bbW^ is around 0.16 fb only (Mh± fti 100 

GeV being a value well below the threshold at uit + mb). The single top resonant 
diagrams via t bW^ (the square of diagram 6 plus its interference with graphs 2 
and 9 in Fig. la) and via t bH^ (the square of the sum of diagrams 1, 4 and 5 
plus their interferences with graph 2 in Fig. la) produce the rates 0.10 and 0.46 fb, 
respectively. At the same point of CM energy the tt contribution is ~ 2.47 fb. As 
^/s increases, the top-antitop signal clearly gets larger, eventually being the dominant 
part of the total rate: e.g., ^ 72 fb at 360 GeV. The Higgs resonances h and A 
yield 0.19 and 0.52 fb at i/i = 360 GeV, respectively, whereas the single-top ones 
produce 0.18 and 0.65 fb, in correspondence of t bW^ and t —>■ bH"^ diagrams. 
The trend remains rather similar for Ma = 100 GeV (and tan/3 = 1.5), whereas when 
Ma = 140 GeV (that is Mh± ~ 161 GeV) phase space suppression effects become 
dominant. On the one hand, the off-shell tt production is well below detection level up 
to ^/s ~ 346 GeV. On the other hand, the h{A) bb resonant graphs are much smaller 
(their BRs remain constant but the sum Mh(A) + Mh± + Mw± approaches (exceeds) the 
CM energy). Therefore, the bulk of the total cross section is due to the signals, 
H"^ hW^ — i> bbW^ being around 2.19 fb and tb ^ bbW^ approximately 3.45 

fb (as Mfj± mt + mh), again at ^/s = 340 GeV. The phase space effect is very drastic 
on top decays, so that the difference between the e'^e" tt curve and that produced 
by all the diagrams is largest in this case (i.e.. Ma = 140 GeV). In fact, an increasing 
value of Mu± means a suppression of the BR[t — > bH"^] and vice versa an enhancement 
of that for if^ —>■ tb bbW^ (also note that for Mh± ^ 161 GeV one gets closer to a 
local maximum of the —>■ hW^ —>■ bbW^ decay fraction |[T7[|). 

For tan/3 = 30. the two H^H~ signals are generally negligible, because of the sup- 
pression coming from their BRs (see Ref. Jl^), the dominant background effects being 
due to the h ^ bb and A ^ bb resonant graphs. At ^/s = 340 GeV and for Ma = 60 
GeV, one gets that hW^ bbW^ and tb bbW^ are at the level of 

(9(10^^) and (9(10~^) fb, respectively, whereas the h{A) —>■ bb resonant graphs yield 
a contribution of 0.39(0.40) fb. However, as Ma increases the h,A ^ bb resonances 
get smaller (they both are around 0.015 fb at Ma = 100 GeV and both completely 
negligible at Ma = 140 GeV, again for y/s = 340 GeV). In fact, although the BRs 
for h and A bosons into bb pairs remain roughly constant for Ma ~ 140 GeV also 
for tan/3 = 30., it should be noticed that M^ grows with Ma more at large than 
at small tan/3's. For example, at Ma = 100(140) GeV and tan/3 = 30., one finds 
Mh + M\Y± + Mh± ~ 308(336) GeV. For such reasons then, the differences between 
e+e" -^tt^ bbW+W- and e+e" bbW+W' rates are small for Ma ^ 100 GeV and 
large tan /3's. 

A second technique of measuring mt at threshold is the so-called end-point method 



in the spectrum of the boson produced in the top decay chain |T^. Whereas for events 
of the type e^e" bbW~^W~ this is necessarily a boson [|, |5|, according to our 
detection strategy such a particle is in events of the type (@)-(0) a if ^ scalar, tagged 
in the rur channel. In general, if one assumes an NWA for the top in e~^e~ — *■ ti (see 
eq. (1) in Ref. then the minimum and maximum values of pn (and pw) are fixed 
by the decay kinematics (see, e.g., p2[)- This is clearly spoiled by a finite value of 
(which introduces a Breit-Wigner dependence in the top propagator), whose effect 
consists in a smearing of the edges of the momentum distribution. The inclusion of all 
the other non-tt diagrams also mimic a similar effect. The pn spectrum in presence of 
both finite top width and all graphs of Figs, la-b is given in the two upper windows 
of Figs. 3a-b (e.g., for = 60 and 100 GeV). For reference, we also give the same 
distribution in the case of the boson (from e'^e" bbW^H^ events), spectrum 
which could well be plotted by means of the momentum of the jet-jet pair for which 
Mjj ^ M]y± . From their comparison with the curves in NWA (lower windows) one can 
easily deduce the difficulty in recognising any end-point in the distributions when all 
effects are properly included, those more visible (the long tails at large momentum for 
e~^e~ —>■ bbW^H"^) being clearly due to the irreducible background. It is also interesting 
to notice the single bin filled on the right hand side of the two upper plots in Fig. 3a, 
for the case of the pn spectrum in reaction (0). This is due to the e+e~ H^H~ 
diagrams (graphs 9 and 12 in Fig. la), for which the momentum of the if^'s is fixed 



via the relation ph = ys/A — Mjj±. The rate of this bin is regulated by the interplay 
between the cross section for e~*'e^ H^H^ and the charged Higgs BRs into hW^ 
and tb pairs, so that the bin practically disappears for tan (3 = 30., for which the latter 
are both heavily suppressed (when Ma ~ 100 GeV). 



yi = 500 GeV 

Far above the 2mt threshold (in our case, y/s = 500 GeV) one can determine the top 
mass from the two three-jet invariant mass distributions that can be reconstructed 
from the bbHrUr final state |T9[. In fact, in bbW^H^ and bbH^H~ events in which one 
final state boson decays hadronically (let us say the positive charged one B~^, where 
= W^, H^) and the other leptonically, there are two possible three-jet combinations 
among the jj pair reconstructing the boson mass (i.e., those for which Mjj ^ Mi4/(i^)±) 
and the two 6's: the 'right' one {B~^b), peaking at rrit, and the 'wrong' one {B~^b), 
producing a rather broad and fiat spectrum (we assume that the charge of the 6-jets 
is not recognised). They are both plotted in Figs. 4a and b, for both reactions (0)- 
d^). For comparison, in case of bbW^H^ events we also give the similar distribution 
which can be obtained by combining the tagged with the two bottom quarks. Note 
that by '6-quarks' we mean here the two jets that do not peak at M^(j:^)±Q. Note 
the logarithmic scale of the figures, as well as the rather thin binning (2 GeV) of the 
distributions. In other terms, although in this case the 'wrong' distributions are well 
below the 'right' ones, it must be remembered that for resolutions worse than that 
considered here the peaks would be cut down, whereas the fiat spectra would remain 
the same. Therefore, the possibility of clearly identifying the top resonance drastically 

^^The tagging procedure can be simplified if /i-vertex tagging techniques are exploited on the two 
&-quarks, so that the requirement Mjj « M^y(^)± needs not to be applied. 



depends on the angle and energy resolution of each single jet. 

For reference, the cross sections at ^/s = 500 GeV for process @ are (no charge 
conjugation included): 105(114) fb for tan/3 = 1.5 and 61(64) fb for tan/3 = 30. when 
Ma = 60 GeV; 61(76) fb for tan/3 = 1.5 and 35(37) fb for tan/3 = 30. when Ma = 100 
GeV. At the same energy, in case of reaction (|^) one obtains: 24(24) fb for tan/3 = 1.5 
and 8.5(8.6) fb for tan (3 = 30. when Ma = 60 GeV; 6.4(6.5) fb for tan/3 = 1.5 and 
2.3(2.4) fb for tan/3 = 30. when Ma = 100 GeV. The numbers in brackets refer to the 
rates obtained when all diagrams are calculated whereas the others to the case of top- 
antitop production and decay only (in both cases F = Ft). From these numbers, one 
realises that in case of e"^e^ hhH^H~ events the effects of the irreducible background 
are negligible, whereas for the process e'^e^ — bbW^H^ they can amount to several 
percent (also at y/s = 500 GeV). 

3.3 Higgs signals 

Let us now turn to the case of Higgs resonances. We have chosen only a few represen- 
tative cases, where Higgs peaks can be visible in the decay spectra and show peculiar 
features. First of all, it is instructive to plot the invariant mass of the bb system, when 
all diagrams are computed. This is done in Fig. 5a for the combination Ma = 60 GeV 
and tan/3 = 1.5, at y/s = 350 GeV, for both processes (|1)-(|1), and in Fig. 5b for the 
case Ma = 140 GeV and tan/3 = 1.5, at = 500 GeV, for reaction (0) only. If 
the resolution in invariant mass of the bb system is 5 GeV, at lower energy only the 
spectrum from process (^ shows a resonance (see Fig. 5a, continuous histogram in the 
main plot). However, this peak actually hides a superposition of two M^^i resonances, 
due to decaying h and A bosons (in order of increasing mass), as it can be appreciated 
in the central insert, where a binning of 1 GeV is used. If such higher resolution can 
be achieved, one can possibly disentangle the h ^ bb peak from reaction (|^) (lower- 
right insert) whereas the H bb one from process (^ remains overwhelmed by the 
background (its position is located by the arrow in the upper-right insert). At larger 
energies, phase space effects do not suppress any longer the production of the heaviest 
neutral Higgs scalar. In fact, at ^/s = 500 GeV, for Ma = 140 GeV and tan /3 = 1.5, the 
resonances of all three neutral Higgses of the AiSSAi can be recognised (see Fig. 5b). 
In general, note that, although such plots might be somewhat superfluous if one wants 
to optimise the cuts in order to enhance e+e~ — >• ti and reduce the background (in fact, 
by the time of the advent of an NLC the values of Mniggs might already be known), in 
contrast their knowledge will be essential if one would like to attempt studies of the 
gauge structure of the AiSSJ^, as mentioned in footnote 6. Like for the case of top 
signals, we again treat the case of charged Higgs signals at y/s = 350 and 500 GeV 
separately. 

= 500 GeV 

We continue this Subsection by discussing the potential of a y/s = 500 GeV NLC in 
detecting charged Higgses produced in e'^e~ H^H^ scatterings, in which one of the 
scalars decays into hW^ — >■ bbW^ and/or tb —* bbW^. In Figs. 6a-b we display the 
differential distributions in the invariant mass of the system bbW^ (i.e., the hadronic 



four-jet subset, as jj), at = 500 GeV and for Ma = 140, 180 and 220 GeV 

(both values of tan/5 are considered). For all pseudoscalar Higgs masses we give the 
results as produced by the full set of Feynman diagrams, whereas in the case Ma = 140 
GeV we also plot the two resonant contributions separately (upper two windows). 
Note that, on the one hand, we consider A masses larger than 140 GeV (i.e., Mh± ~ 161 
GeV) to allow for large values of the two Higgs BRs into hW^ and tb (either or both), on 
the other hand, we restrict ourself to the case Ma ~ 220 GeV because of the maximum 
CM energy available. It is interesting to notice how the spectra produced by the non- 
H^H~ diagrams steeply increase at M^g^i/i ^ 170 — 180 GeV (roughly, the top rest 
mass: compare the lower window in Fig. 6a against the two upper ones). This is 
somewhat fortunate, as for Ma = 140 GeV the charged Higgs peak lies right a few 
GeV below that. Furthermore, for this value of Ma and for tan/? = 1.5, the Higgs 
production and decay rates are similar in the two channels (compare the continuous 
curves in the two upper frames of Fig. 6a). This could allow one (after implementing 
an appropriate selection requirement: e.g., M^j f» Mh) to separate the hW decay from 
the th one. In fact, for larger values of Mh± the Higgs signals come in practise from 
top-bottom decays only (see Fig. 6b). For tan/5 = 30. this is also true at Ma ~ 140 



GeV ||T^, as it can be appreciated by comparing the dashed curves in the two upper 
sections of Fig. 6a. 

Since the rates given in Figs. 6a-b do not include the BR into rvr pairs of one of 
the H"^ bosons and since this becomes very small for Ma ~ 140 GeV at small tan/3's 
1l7| , one has to face the following scenario (assuming the nominal integrated luminosity 



/ Cdt = 10 fb~^). If tan/? is large, charged Higgs signals can be produced at large rate 
at an NLC with y/s = 500 GeV and certainly detected in the H —>■ tb ^ bbW^ channel 
for Mh± ~ 161 GeV. In contrast, for small tan/?'s, this is true only for Ma ~ 140 — 150 
GeV, region where the 'narrow' contribution from the hW channel compensates the 
small below threshold branching ratio into tb pairs. 

= 350 GeV 

Concerning smaller values of CM energy, that is y/s = 350 GeV, it must be said that 
phase space suppression largely counterbalances the increase of the production cross 
section due to the s-channel dynamics of e~^e~ H^H~ events. In the very end 
the spectrum of charged Higgs masses which can be covered by the two mentioned 
decay channels at that energy is rather narrow, about 20 GeV, and only if tan/5 = 
1.5. In fact, for such a value, when Mh± is below 130 GeV or so, the total BR into 
the two Higgs channels is smaller than 0.1%. When instead Mh± ~ 150 GeV the 
production cross section fall below detection level |T^. If tan/3 = 30. such a window 
disappears completely. For these reason then, we decided not to focus our attention 
on the threshold energy stage of the NLC, in the case of charged Higgs searches via 
e+e" ^ 



3.4 Mass dependence 



So far, we have concentrated in our analysis on the case of selected values of Ma 
(or, equivalently, of Mh±)- However, as the effects due to the finite widths and to 



the irreducible backgrounds are equally present over all the available charged Higgs 
mass range, in this closing Subsection we intend to generalise our results by presenting 
various rates for processes (|^)-(|^) as a function of Mij± (again, at fixed tan/3 = 1.5 
and 30.). We do so, e.g., at the CM energy value of ^/i = 500 GeV, as we have noticed 
that a second stage NLC allows one to study contemporaneously both e+e~ — >• ti and 
e^e~ —>■ H^H events. 

Fig. 7 shows the cross sections for the same three processes as defined in Tab. I-II, for 
both reactions but now at ^/s = 500 GeV and over the continuous mass range 

60 GeV ^Ma^ 240 GeV (i.e., 100 GeV ^ Mh± ^ 252 GeV). A comparison between the 
rates obtained in NWA (i.e., continuous lines) and those from the production and decay 
diagram only (i.e., dashed lines) illustrates that finite width effects clearly becomes 
larger as the difference AM = rrit — Mh± — ruh gets smaller, the latter being positive. If 
negative, the rates of the NWA are identically zero. For Mfj± well below rrit the effects 
are negligible (also compare to Figs. 2a-b at y/s = 350 GeV). Furthermore, as the top 
width does not vary drastically, from approximately 1.72(1.78) GeV at = 60 GeV 
and tan (3 = 1.5(30.) down to around 1.40 GeV at Ma = 150 GeV, where only the SM 
decay channel t —>■ hW^ survives (see Fig. 2b of Ref. 0), width effects are (above the 
t — >■ hH^ threshold) rather independent from Mu±. The total rate of process (H) is, in 
comparison to that of e^e~ ^ ti ^ bbW^H^ events, much larger than the cross section 
of reaction (|^) is with respect to the corresponding e~^e~ tt bbH^H~ rates. As 
already outlined in Section 3.1, this is due to the much richer variety of non-tt events 
that are active in the first case (especially the charged Higgs resonances: notice, e.g., 
the onset of the top-bottom channel at large M/^±'s and tan/5 = 30.). 

In Fig. 8 we have plotted the cross sections of the two possible Higgs mechanisms en- 
tering in e+e" — > bbW'^H'^ events: i.e., via if^ — ^ tb and iJ^ hW^ decays (actually, 
the second channel is only considered at small tan /3's, because of the vanishing BR at 
larger values fl^])- For both channels, the figure shows rates as obtained by multiplying 
the total e^e~ — *■ H^H^ production rates times the branching ratios BR(if^ tb) and 
BR(if^ — >■ hW^) (i.e., assuming on-shell production and neglecting the Higgs widths) 
and those produced by the subamplitudes involving the two mentioned resonances in 
our ME, along with the yield of all e^e~ —>■ bbW^H^ diagrams. Whereas the effects of 
the finite width of the i7^-bosons are generally small in both Higgs channels (possibly 
apart from the hW^ channel for Mh± ~ mt at tan/3 = 1.5), the irreducible 

background is overwhelming the signals up to Mh± ~ 150 GeV (that is, as long as 
the t — >■ bH"^ channel is well open), the rates of the two remaining comparable up to 
charged Higgs masses around 200 GeV. Above Mh± ^ 150 GeV, it is the H"^ — >• tb 
component of the e~^e~ — > bbW^H"^ cross sections which accounts for the total rate, 
this being modulated at large values of Mh± by a phase space suppression. Since the 
largest part of the cross section of non-Higgs diagrams in e~^e~ —>■ bbW^H^ events is 
due to tt events, the implementation of a cut like Mjjj ^ rrit (on the three-jet invariant 
mass distributions that can be reconstructed from the fcfejjrz/,- signature) should allow 
one to largely increase the Higgs signal-to-background ratio. However, notice that in 
Section 3.3 no such cut was enforced, the Higgs peaks still being visible in many in- 
stances. This makes clear that concrete chances of Higgs detection exist for charged 
scalars of the AiSSAi produced via e'^e~ H^H~ interactions and decaying via 
either hW^ or ^ tb ^/s = 500 GeV e+e" Linear CoUiders. 



4. Summary and conclusions 



In this paper we have studied the processes 

e+e- ^ bbW^H"^ and e+e" ^ bbH+H-, 

by computing their complete matrix elements at tree-level and integrating these over 
the corresponding phase spaces to produce total and differential cross sections relevant 
to phenomenological studies at the NLC. We focused our attention on the CM energies 
of ^/s = 350 GeV (top-antitop threshold, first stage NLC) and 500 GeV (final stage 
NLC). 

These processes are experimentally interesting since they both include among the 
various contributions the production of top-antitop pairs eventually decaying via the 
channels t — > bW^ and t — > feif^. Furthermore, in the case of bbW^H"^ final states one 
can also study the charged Higgs production process e~^e~ H^H~ followed by the 
decay of one of the if^'s via hW and/or tb pairs, with h bb and t bW. 

Numerical analyses have been carried out in order to quantify, first, the influence 
of finite top width and irreducible background effects on the integrated and differen- 
tial rates as obtained from e+e~ — > tt events in NWA and, second, to estabhsh the 
detectability of the mentioned Higgs processes, as a function of the values assumed by 
the fundamental parameters of the AiSSAi. 

For certain combinations of Ma and tan/3, our results indicate that, on the one 
hand, finite width and background effects must be included in the simulations aiming 
to measure precisely the top parameters and, on the other hand, charged Higgs signal 
in the mass range 140 GeV ~ Mh± ~ V^/'^ can be promptly detected at nominal 
luminosity, especially if ^/s and tan/3 are large. 

The numerical codes we have used in our analysis can be easily folded with the 
experimental Monte Carlo simulation programs (including beam and detector effects) 
used by the NLC Working Groups, in order to quantify correctly irreducible background 
effects in top studies as well as to optimise the selection strategy of charged Higgs 
signals. In this respect, note that we have always assumed one of the two bosons in 
the final state to decay into pairs, which is clearly restrictive in the case of Higgs 
searches for large values of Ma and small tan/3's. Indeed, dedicated experimental 
simulations might well assess that other H^H" decay signatures are also relevant, so 
that we make our matrix elements available to the public for further studies. 

Finally, we also have pointed out the importance of events of the type discussed 
here in order to study the gauge structure of the Higgs sector of the M.SSM., as they 
involve many of the tree-level Higgs vertices of the theory. 
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Table Captions 



Tab. I Cross sections in femtobarns for e^e~ — > tt (NWA, F — > 0), for e^e~ tt ^ 
bbW^H^ (production and decay diagram only, F = Ft) and for e+e~ — > bbW^H^ 
(all diagrams at trcc-lcvcl, F = F^), at ^/s = 350 GcV, for mt = 17 A GeV, and 
Ma = 60(140)[220] GeV (i.e., Mh± « 100(161)[234]). The first(second) row 
reports rates for tan/3 = 1.5(30.). 

Tab. II Cross sections in femtobarns for e'^e~ — >■ tt (NWA, F — >■ 0), for e'^e~ ^ ti ^ 
bbH'^H~ (production and decay diagram only, F = F^) and for e~^e~ bbH^H~ 
(all diagrams at tree-level, F = Ff), at ^/s = 350 GeV, for rrit = 174 GeV, and 
Ma = 60(140)[220] GeV (i.e., Mh± ~ 100(161)[234]). The first(second) row 
reports rates for tan/3 = 1.5(30.). 



Figure Captions 



Fig. 1 The AiSSAi Feynman topologies contributing at lowest order to process (|^) 
(a) and (§) (b). The label 'h(h')' represents any of the neutral Higgses of the 
theory, H,h and A, whereas 'H' refers to the charged scalar, H^. 'A', 'Z' and 
'W represent the gauge bosons 7, Z and W^. Note that for process (|) we plot 
the diagrams for the charge combination W~H~^ only, those corresponding to the 
case W~^H~ being trivially deducible from the one given here. The figure has 
been produced with the help of the package MadGraph . For reference, in (c) 



we present the branching ratios of t-quarks and if ^-bosons as functions of the 
charged Higgs mass in the range 60 GeV ~ Mh± ~ 240 GeV, for tan/? = 1.5 and 
30. 

Fig. 2 Cross section in femtobarns for e^e~ bbW^H"^ (a) and e^e~ — > bbH~^H~ 
(b) events around the top-antitop threshold, for three different values of Ma- 
Continuous lines: tan/3 = 1.5 and production and decay diagram only. Dashed 
lines: tan/3 = 1.5 and all diagrams. Dotted lines: tan/? = 30. and production 
and decay diagram only. Dash-dotted lines: tan (3 = 30. and all diagrams. In all 
cases r = has been used. Note the superposition of curves in some instances. 

Fig. 3 Differential distributions in the momentum of the and bosons, in 
e^e~ —>■ bbW^H^ and e+e~ —>■ bbH^H^ events, for two values of Ma, with 
tan/5 = 1.5 (a) and 30. (b). The two following cases are considered: NWA (lower 
windows) and all diagrams (upper windows). The CM energy is ^/s = 350 GeV. 
Continuous lines: momentum in e~^e~ —>■ bbW^H"^ events. Dashed lines: 
momentum in e+e~ — > bbW^H^ events. Dotted lines: momentum in 
e'^e~ — >■ bbH~^H~ events. 

Fig. 4 Differential distributions in the invariant mass of the 'right' and 'wrong' three- 
jet combinations (see the text), in e~^e~ bbW^H"^ and e~^e~ bbH^H~ 
events, for Ma = 60 (upper window) and 100 (lower window) GeV, with tan (3 = 
1.5 (a) and 30. (b). All diagrams have been here considered. The CM energy 
is a/s = 500 GeV. Continuous lines: Mhw± mass in e+e~ — >• bbW'^H'^ events. 
Dashed lines: Mi,h± mass in e+e" — >• bbW^H^ events. Dotted lines: M^h^ mass 
in e'^e~ — > bbH~^H~ events. 

Fig. 5 Differential distributions in invariant mass of the bb pair, (a) In e^e~ — >■ 
bbW^H^ and e+e" ^ bbH+R- events, for tan/? = 1.5 GeV, Ma = 60 GeV 
and ^/s = 350 GeV. Continuous lines (main window): all e"'"e~ bbW^H^ 
diagrams. Dashed lines (main window): all e+e~ bbH^H^ diagrams. The two 
small windows on the right show the regions around the H (in e+e~ — > hbW^H^, 
upper one) and h (in e^e~ — > bbH~^H~, lower one) resonances enlarged and 
plotted with high resolution, together with the tt contribution (F = Ff, shaded). 
The small central window does the same for the h and A 'resolved' resonances 
in the case of e~^e~ — > bbW^H^ events, (b) In e"'"e~ bbW^H^ events, for 
tan/3 = 1.5 GeV, Ma = 140 GeV and ^/s = 500 GeV. Dashed line: all diagrams. 
Continuous line: ti contribution (F = F^, shaded). 



Fig. 6 Differential distributions in invariant mass of the bbW^ system, in e+e — > 
bbW^H^ events, for tan/3 = 1.5 and 30., Ma = 140 (a) and 180, 220 GeV (b), 
at ^ = 500 GeV. In the case = 140 GeV the resonant hW^ (upper 

plot) and if^ — > tb (middle plot) contributions are shown separately. Continuous 
lines: tan/3 = 1.5. Dashed lines: tan/3 = 30. 

Fig. 7 Cross section in femtobarns for e^e~ bbW^H^ (upper three curves) and 
e+e" — >• bbH~^H~ (lower three curves) events at y/s = 500 GeV as a function 
of Mh±, for tan/3 = 1.5 (upper plot) and tan/3 = 30. (lower plot). Continuous 
hues: e+e" tt (NWA, T ^ 0). Dashed hues: e+e" ^ ti ^ bbW^H'P and 
e+e" — > — > bbH^H~ (production and decay diagrams only, F = F^). Dotted 
hues: e+e~ bbW^H^ and e+e" bbH^H~ (all diagrams at tree-level, F = 
Tt). 

Fig. 8 Cross section in femtobarns for e"*"e~ bbW^H^ events at y/s = 500 GeV 
as a function of Mh±, for tan/3 = 1.5 (upper plot) and tan/3 = 30. (lower plot). 
Continuous lines: cr(e+e" — > H^H~) x BR(if^ ^ tb ^ bbW^). Dashed lines: 
e+e- ^ H+H~ tbH"^ bbW'^H'^. Dotted lines: a{e+e- H+H-) x 
BR{H^ hW^ bbW^). Double-dash lines: e+e" ^ H+H' hW^H^ 
bbW^H^. Dash-dotted hues: e+e" — > bbW^H^ (all diagrams at tree-level, F = 



mt (GeV) 


a{e+e- X) (fb) 




tt 


tt bbW^m 




174 


34.64(2.21) [1.32 x 10"^^] 
20.04(1.30) [6.37 X lO"^^] 


30.55(1. 73)[2.57x 10"^] 
17.57(1.01)[1.24 X 10-^] 


40.50(10.32) [9.48 x 10"^] 
18.89(1.05)[9.33 x 10"^] 


^ = 350 GeV 



Tab. I 



mt (GeV) 


cr(e+e- ^ X) (fb) 




tt 


tt bbH+H- 


bbH+H- 


174 


7.87(2.17 X 10-2)[o.OO] 
2.84(7.50 X 10-=^) [0.00] 


6.98(1.46 X 10-2)[0.00] 
2.49(5.04 X 10-3) [0.00] 


7.10(1.46 X 10-2)[0.00] 
2.52(5.05 X 10-3) [0.00] 


= 350 GeV 
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Diagrams by MadGraph 
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